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Abstract 

o 

"p^ I We construct hadronic axion models in the framework of the anomaly-mediated 

supersymmetry breaking scenario. If the Peccei-Quinn symmetry breaking is related 
^ , to the supersymmetry breaking, mass spectrum of the minimal anomaly-mediated 

Qj I scenario is modified, which may solve the negative slepton mass problem in the min- 

>-Ch ■ imal anomaly-mediated model. We find several classes of phenomenologically viable 

models of axion within the framework of the anomaly mediation and, in particu- 
^ \ lar, we point out a new mechanism of stabilizing the axion potential. In this class 

;_j ' of models, the Peccei-Quinn scale is related to the messenger scale. We also study 

phenomenological aspects of this class of models. We will see that, in some case, 
the lightest particle among the superpartners of the standard-model particles is stau 
while the lightest superparticle becomes the axino, the superpartner of the axion. 
With such a unique mass spectrum, conventional studies of the collider physics and 
cosmology for supersymmetric models should be altered. 



1 Introduction 



In modern particle physics, symmetries play extremely significant roles in solving prob- 
lems related to various fine tunings. Among them, in this paper, we consider two serious 
problems, that is, the strong CP problem and the gauge hierarchy problem. 

Taking account of the instanton effect, the B-parameter in QCD, which is given byp^ 

with being the gauge coupling constant for SU{?>)c and GJJ^ the field strength of the 
gluon, is constrained to be smaller than 10^^ otherwise, the electric dipole moment of 
the neutron becomes too large to be consistent with the experimental constraint. Such a 
smallness is, however, unnatural in the standard model, since is just a parameter in the 
Lagrangian. One smart solution to this problem is to introduce the Peccei-Quinn (PQ) 
symmetry |]^, ^ ^, which is a spontaneously broken global abelian symmetry which is 
anomalous under QCD. With such a symmetry, which we call f/(l)pQ, Nambu-Goldstone 
boson (called "axion" a) shows up and couples to the gluon as 

where /„ is the decay constant of the axion. In this case, the parameter is promoted to 
a dynamical variable. Minimizing the potential of the axion, 0efT = + a//a — >^ 0, and 
hence there is no strong CP problem in this case. Since the PQ symmetry solves the strong 
CP problem very beautifully, it is desirable to implement the PQ symmetry in models of 
new physics beyond the standard model. 

Another serious fine-tuning in the standard model is related to the stability of the 
electroweak scale against radiative corrections; in the standard model, radiative corrections 
to the Higgs mass parameter quadratically diverge and hence the electroweak scale is 
expected to be as large as the cutoff scale. Taking the cutoff at the gravitational scale, we 
need more than 30 orders of magnitude of fine tuning. 

Once the supersymmetry (SUSY) is introduced, quadratic divergences cancel between 
bosonic and fermionic diagrams. Thus, SUSY is currently regarded as one of the most 
promising candidate of the solution to the naturalness problem, and we consider axion 
model in supersymmetric framework.P^ Among various models, in this paper, we consider 



the anomaly- mediated SUSY breaking model [0]. Interestingly, the anomaly-mediated 
model may solve the SUSY CP and FCNC problems without fine tuning, since in this 
model, all SUSY breaking parameters are well controlled by the super- Weyl (SW) anomaly. 
In addition, this model may provide solutions to cosmological difficulties such as gravitino 
problem and cosmological moduli problem |^. Thus, we will study how axion models 

The Q parameter here is defined in the bases where all the quark masses are real. 

For early attempts to construct axion models in the supersymmetric models, see Refs. B 
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Table 1: Particle content of the model. The index i is the flavor index which runs from 1 

to A^5. 



can be constructed in the framework of the anomaly-mediated models. In particular, in 
some class of the axion model, we will emphasize that the axion multiplet changes the 
prediction of the minimal anomaly-mediated model and that the negative slepton mass 
problem may be solved. We will show that, in this case, the resulting superparticle mass 
spectrum is similar to that of the deflected anomaly mediation |[TD|. 



Organization of this paper is as follows. In Section 0, we introduce several models of 
axion in the framework of the anomaly-mediated SUSY breaking. Then, in Section ^, we 
study phenomenology of models introduced in Section ^ We summarize our results in 
Section ||. 



2 Model 

Let us discuss the framework of our model. Although our mechanism works in a large 
class of models with various particle content, we consider a model with A^^s-pairs of 5 and 
5 representations in the SU{5) as an example. The particle content of our model is given 
in Table where the representations for the standard-model gauge groups as well as the 
charge for the U{1)pq are also shown. In our model, the lowest component of X acquires 
a vacuum expectation value (VEV) and it breaks the f/(l)pQ symmetry. Thus, we call X 
as the axion multiplet. 

In addition to the fields in the observable sector, we also introduce hidden-sector field 
which is responsible for the SUSY breaking. We denote the hidden-sector field as z. In our 
model, we adopt the sequestered structure 



K = -3log[az\z)+ax\X)], (2.1) 

which may arise, for example, when the hidden and observable sectors are geometrically 
separated |ll2[ . 



Once the sequestered structure is assumed, phenomenology in the observable sector is 
insensitive to how the SUSY is broken in the hidden sector since the effect of the SUSY 



The function ^ also depends on the observable sector fields such as quark, lepton, and Higgs multiplets. 
Such fields are irrelevant in studying the axion potential, and we omit these fields in the expressions. 
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breaking is mediated to the observable sector only by the SW anomaly. Thus, we only 
assume that the SUSY is somehow broken in the hidden sector. In this framework, SUSY 
breakings in the observable sector are only from couplings to the compensator field $, 
whose VEV is given by 



$ = 1 + 



[2.2) 



Here, F$ is an auxiliary field in the gravitational multiplet, and its VEV is given by, 
assuming vanishing cosmological constant. 



1 IdK ^ 



(2.3) 



where ~ 2.4 x 10 GeV is the reduced Planck scale. With this compensator field, 
relevant part of the supergravity Lagrangian has the following form [jl3 

C - j d^e^^^e-^^'^ + j d^e^^W + h.c 



(2.4) 



Soft SUSY breaking parameters in the observable sector are obtained by expanding the 
above Lagrangian in the background with non- vanishing F<^. 

Now, let us consider the potential of X and study how the PQ symmetry can be 
broken.p^ As we will see below, there are a couple of different approaches, each of which 
have different consequences on the mass spectrum of the squarks and sleptons. Thus, we 
will study each of them separately. 

Importantly, because of the [/(1)pq symmetry, the function ^ in the Kahler potential 
depends on the combination X'^X. At the tree level, expanding ^ around \X\ = 0, we 
express 



e,.ee(Xt, X) = Ml c, (^] ' = CiXtX + ^(XtX)2 + 



p=l 



(2.5) 



where Cp are constants. Notice that, as we will see below, radiative correction to the 
function ^ becomes significant in some case. In addition, the axion multiplet X is coupled 
to other fields in the superpotential as 



W^X55 = AqXQQ + \lXLL. 



(2.6) 



Due to this superpotential, the axion multiplet X couples to the field strength of the gluon 
multiplet, and the imaginary part of the lowest component of X becomes the axion. 



We use the same expression for the chiral superfield and its bosonic (lowest) component. 
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At the tree level, and with the simplest Kahler potential of ^ = ciX'^X, VEV of X is 
undetermined. Indeed, after the rescaling 

X = ^X, Q = ^Q, (2.7) 

all the compensators in Eq. (|2.4| ) are absorbed in the "hatted" superfields. (Hereafter, 
"hatted" superfields denote superfields after the rescaling with the compensator like 
Eq. ( p.7|) .) In this case, no potential is generated for X. Taking account of the higher 
dimensional operators and/or radiative corrections, however, non-trivial potential for X is 
generated. In the following subsections, we study several cases. 



2.1 Model 1 



The first approach is to use the radiative correction to the wave function renormalization 
factor of X and the higher dimensional terms in the Kahler potential. Taking account of 
the wave function renormalization, at the loop level, we expect the function ^ to become 



ax\ X) = c,zx{x\ x)x^x + ^{x^xy + 



(2. 



The X-dependence of the wave function renormalization factor Zx is from the scale de- 
pendence; below the scale of the VEV of X, all the particles coupled to X decouple and 
Zx does not run. Thus, we obtain 



zx{x\x) = y: 
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1, x^xV 



A2 ] 



+ -Z^(A) -log 



xtxV 



+ 



(2.9) 



where A is an arbitrary scale, and the "prime" represents derivative with respect to log fi. 

Taking account of these effects, the supergravity Lagrangian given in ( |2.4|) becomes, 
after the rescaling (p.7|). 



C 



d^ez 
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2Zx $^<I> 



X^X IZ'l / x^xX 



'X 

^Zx 



+ ■ ■ 



X^X 



C2 {X^Xf 
M2 <I)t$ 



+ 



(2.10) 



where we performed an X-independent rescaling of X, and the potential for X is given by 



V{X\X) 



m 



d{\ogfif 



\F^\^X^X 



J fj.=XX 



_C2_ 
M2 



\F^\\X^Xf + 



XtX-^|F,p(XtX)^ + 



(2.11) 
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With the superpotential given in Eq. ( |2.6| ), the mass for the X field is given by, at the 
one-loop level, 



m 



X 



(Xt,X) 



X 



0?(log/i)2 



J fi=XX 



167r2 



iQgi + -gt]Xt3 + [Q92 + '^9t)K 



-iV5(3A^ + 2Ai)^-2(3A^ + 2Ai) 



16gl - 5(5A^5 + 2)A^ 



(2.12) 



where, in the last equahty, we used the approximations 3> g2 and gi, and A = Aq ~ A^. 
Hereafter, we adopt these approximations for simplicity. Our results are qualitatively 
unchanged even if we use the exact formula. 

For the spontaneous breaking of the PQ symmetry, should become negative. This 
happens when the gauge coupling constants are larger than the Yukawa couplings \q 
and Xl- Such a situation can be naturally realized by assuming relatively small values 
of the coupling constants Xq and A^,. Of course, for a realistic model of SUSY and PQ 
symmetry breakings, potential of X should be somehow stabilized. Pomarol and Rattazzi 
|10| proposed to use the inverted hierarchy mechanism [|ll| to stabilize the potential; if 
m\ changes its sign at some scale with drrij^ /d log fi > 0, the potential of X has stable 
minimum at the scale = 0. Such a scenario is possible when there exists asymptotically- 
free gauge interaction with relatively large gauge coupling constant. We found, however, 
that the SU{3)c gauge interaction cannot play this role since its gauge coupling constant 
is not large enough to realize the inverted hierarchy mechanism. We checked that this 
is always the case, by solving renormalization group equations numerically. Pomarol and 
Rattazzi suggested to gauge (some part of) the fiavor symmetry to introduce a new gauge 
interaction for the stabilization. Unfortunately, it does not work for the axion model since, 
once Q and Q (and/or L and L) have non-trivial gauge quantum numbers under extra 
gauge interaction, the axion potential is modified and the strong CP problem cannot be 
solved. 

In our model, we do not rely on the inverted hierarchy mechanism, since there is another 
simple way to stabilize the potential. As suggested in Eq. (p.ll|) , the higher dimensional 
terms in the Kahler potential naturally exist which result in quartic and higher order terms 
in the scalar potential. Assuming the negativity of C2, the potential given in Eq. ( p.ll| ) has 
a minimum, and the VEV of X is given by 



J a 



m 



X 



2cc 



Mt 



gsX j r~ 

47r2V2|c2|Ar5 



(2.13) 



where we assumed g^ ^ A. As one can see, the PQ scale is suppressed by the loop factor, 
and it decreases as the coupling constant A becomes smaller. Even with A ~ O(10~^), /„ 
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is 0(10^^ GeV), and fa becomes 0(10^^ GeV) when A ~ 0(10^'^). Phenomenological and 
cosmological imphcations of this fact will be discussed later. 

Next, let us discuss the mass spectrum of the superparticles in the MSSM sector. Since 
the X field couples to the chiral multiplets which have standard- mo del quantum numbers, 
it affects the soft SUSY breaking parameters in the observable sector at the loop level. 
Indeed, at the scale of the VEV of X, the chiral superfields Q, Q, L and L decouple. 
Then, the running effects from the cutoff scale A to the scale of the VEV of X is via the 
combination log(X'fX/A^) = \og{X^X /A'^^^^). Since the VEV of the highest component 
of X (i.e., the axion multiplet after the rescaling), which is denoted as Fj^, vanishes, we 
obtain F$ dependence from log(X^X/A^$^$). Expanding this logarithm, we obtain the 
SUSY breaking parameter at the messenger scale 



Ma(M„ 



'ai(M„ 



(2.14) 



1 



(47r)2 
-6ai(Mme 

y/(Mmesi 

47r 



2Ci {h-N^)al{M^,,,) 




ai(M^ess) + 3a2(M„ 



16 



a3(M„ 



where the messenger scale is given by 



2C(ai{M^,,,)-N^at{M^ 



A(X), 



(2.15) 



(2.16) 



(2.17) 



and at is top-quark Yukawa coupling, ctj are gauge coupling constants {i runs over the 
MSSM gauge groups.), h = (— y,— 1,3), and C-^ is the second-order Casimir. (For the 
fundamental representations of SU{3)c and SU{2)l, = | and C2 = |, and for f/(l)y, 
C( = with Y being the hypercharge quantum number.) In addition, Nu = (1, 2, 3) for 
qj^'^, tR, and and = for other particles. 

Contrary to the conventional scenario, there is another important particle in our model, 
that is, the axino which is the superpartner of the axion. The axino may become the LSP 
in our model which has significant implications for collider physics and cosmology. Axino 
mass arises from the Lagrangian ( |2.1CI| ); expanding the Lagrangian, the axino mass is given 
by 

1 £Zx (Xt) . ^.2A2 

^'^ = -aUT^IVV^^ ^ T^^si^*' 2.18 
4ct(logyu)^ (A) 167r'* 

where we neglected unimportant phase in the second equality. As one can see, the axino 
mass arises at the two-loop level and is much smaller than the masses of the superparticles 
in the MSSM sector. Thus, in the model 1, the axino becomes the LSP. 
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Interactions of the axino with the observable sector fields depend on the PQ charges of 
the Higgs (and other) fields. The charge assignment is model-dependent, and in particular, 
the charge assignment is related to how the ij,h- and S^-parameters are generated. This 
issue will be discussed later. 



2.2 Model 2 



It is also possible to break the PQ symmetry by higher dimensional operators in the super- 
potential. For this purpose, we introduce extra chiral multiplet Y. With this superfield, 
we write down the following superpotential 



(2.19) 



where px and py are positive integers, and A is some mass parameter which is assumed to 
be of order the Planck scale. Obviously, this superpotential is C/(1)pq invariant assigning 
the PQ charge —px/pv to Y. 

In the background with non- vanishing F$, the scalar potential is given by 



V = 



2 r 



+ {Px+Py - 3) 



j^Px-iyP"' 



+ Py 



\px+PY-3 



XPx yPY-^ 
XPxyPY + h.c. 



Minimizing this potential, we obtain 

Py\{X)\' = px\{Y)\' = C 

where 

2/{px+PY-2) 



2/{px+PY-2) 



^ ^ hx+PY-^ \' 

\Px+Py - 1/ 



^(pY-2)/ipx+PY-2)(px-2)/{px+PY-2) 

Px Py 1 



which is a constant of 0(1). Thus, we obtain 

1(f) I \F^\"-^f^"''^ 



fa 



and 



X 



p-3 
'p-1 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



where p = px + Py- In this case, the PQ scale is determined by the relative size of F$ 
and A. Taking A ~ M*, for example, and F$ ~ 100 TeV, we obtain ~ O(10^° GeV), 



0(10^^ GeV), and 0(10^^ GeV), for p = 4, 5, and 6, respectively. Especially, in the case 
ofp = 4, N^>10. 



With the VEV given in Eq. ( |2.24| ), we obtain 
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(2.25) 
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(2.26) 
(2.27) 



One important difference between this model and the model 1 is the axino mass. In this 
model, axino is a linear combination of the fermionic component in the chiral superfields 
X and Y. The mass matrix for these fermions are, after relevant phase rotation to remove 
unwanted phases. 



axino 



(jipx+PY)/2pi2'PY)/2^{2-px)/2 f P X - 1 ^JVXVY\^ 



Py 



V VPxPy Py - I ) 



\F^\. 



(2.28) 



As one can easily see, both fermions acquire masses as large as \Fq,\, which is of order 
10 — 100 TeV. Thus, in this model, axino mass is much heavier than the masses of the 
superparticles in the MSSM sector and hence the axino cannot be the LSP. 



2.3 Model 3 

In the previous models, the PQ symmetry is unbroken in the supersymmetric limit. In some 
case, however, the PQ symmetry can be spontaneously broken even if the supersymmetry 
is preserved. Let us finally consider such 

The simplest superpotential realizing such a situation is 

W = WxE5 + KY'iXX - Nif!), (2.29) 

where fa here is some constant, and Y', whose U{1)pq charge is 0, is a chiral superfield 
which is singlet under SU{5). Contrary to the previous models where the PQ scale fa is 
somehow related to the SUSY breaking parameter /„ is an input parameter in this 
case. As one can see, the PQ symmetry is broken solving the condition for the vanishing 
F-component of Y': dW/dY' = 0. 

In this model, in fact, the PQ symmetry is broken but the masses of the squarks and 
sleptons are not modified. This can be easily seen by solving the equations of motion. 
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After the rescalings X = $X and so on, we obtain {XX) = f^^^, by solving dW/dY' = 0. 
At this stage, the relative size of X and X is undetermined since it corresponds to a flat 
direction in the supersymmetric limit. The relative size is determined taking account the 
effect of the SUSY breaking. With non- vanishing Y' acquires a VEV as 

k{Y') ~ (2.30) 

which gives equal masses to X and X. Then, the VEVs of these chiral superfields become 
the same and 

(X) = (X) = Ar5/,$, (2.31) 

up to an irrelevant phase. Thus, substituting this VEV into log(X^X/A^$^$), effect of 
the SUSY breaking disappears, and hence the axion multiplet does not modify the masses 
of the superparticles. In this case, extra mechanism should be introduced to solve the 
negative slepton mass problem, and we do not pursue this direction anymore. 



2.4 The fiH- and 5^^- Parameters 

Before closing this section, let us comment on the fin- and 5^-parameters. The /ih- and 
i?^-parameters can be generated by slightly modifying the model by Pomarol and Rattazzi 
10| . The new superpotential we introduce is 



W = XsTxSTX + XsstS'T + XthhTH^H,, (2.32) 

where 5* and T are chiral superfield which are singlet under the SU{5), whose PQ charges 
are +1 and —2, respectively. Notice that we also assigned the PQ charge —1 for the up- 
and down-type Higgses. 

As in the case of the model by Pomarol and Rattazzi, kinetic mixing between X and 
S appears at the one-loop level: 

Cmi. = J d^e<l>^<l>Zsx{fi = \XX\)SX^ + h.c. (2.33) 

In addition, once X acquires a VEV, S and T becomes massive and these fields can be 
integrated out. In particular, by solving the condition dW/dT = 0, we obtain 

^^_W^^ (2.34) 

XsTX X 



and substituting Eq. ( p.34| ) into Eq. ( p.33| ), we obtain 



Crni. = / d'eZsxifi = \XXm)^HuHd + h.c. (2.35) 

^STX J X 
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The model 1 



The model 2 (p = 4 case) 
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Table 2: Mass spectra of the model 1 and 2 in units of GeV. We take tan (3 = 30, M^ess = 
lO^^GeV and lO^^GeV, Nr, = 7 and 13 for the model 1 and 2, respectively. The Wino mass 
M2 is taken to be 600 GeV for the model 1 and 1450 GeV for the model 2. 



In Model 1 where Fj^ = 0, expanding the above Lagrangian, we obtain the fin- and 
parameters to be 
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IJL=\X 



3 Phenomenology 

Here we would like to discuss phenomenological issues of the models we presented. As we 
mentioned earlier, our models have superparticle masses identical to those in the deflected 



anomaly mediation. Some of our results were already discussed in Refs. |10, 14 



3.1 Mass Spectrum 



Let us first consider the model 1. In this model, the soft masses are given in Eqs. ( p.l4| ), 
( p.l5| ) and ( p.l6| ), which are parameterized by the scale of the SUSY breaking the 
messenger scale Mmess and the number of the messenger fields N^. In addition, we have the 
liH and Bfj^ parameters. Here we do not take any particular mechanism to generate them. 
fiH is solely determined so as to reproduce the correct electroweak scale, and B^ is related 



10 



1014 p 

1013 . 




Figure 1: Lower bound on the messenger scale Mmess as a function of the number of 
messengers N^. The hght-shaded region is excluded by tachyonic slepton for tan/3 = 5. 
The lower limit for tan/3 = 30 is also shown as a solid line. The dark-shaded region is 
excluded because the gauge coupling constants blow up below the GUT scale. 



to tan the ratio of the vacuum expectation values of the two Higgs in the MSSM. Thus 
the superparticle mass spectrum is specified by Mmess, and tan/3. 

To obtain a realistic mass spectrum (positive slepton masses etc.), the number of the 
messenger should be large enough. In fact for = 3, the gluino mass vanishes, resulting 
in an unrealistic mass spectrum. For = 4, we found that the sleptons are too light 
to survive the experimental mass bounds, unless the SUSY breaking scale is very high. 
This is disfavored, because the other superparticle masses will exceed 1 TeV, causing the 
fine-tuning problem. Thus we will consider the case A^s > 5. In Fig. |1], we show the lower 
bound on the messenger scale Mmess as a function of N^. The bound comes from the two 
requirements: (i) the positive slepton mass squared (ii) the perturbativity of the gauge 
coupling constants up to the GUT scale. The former requirement is severer for smaller N^, 
while the latter is severer for larger Nr,. One finds that the messenger scale Mmess must be 
larger than ~ 10^ GeV.p^ 



In Figs. and ^ the masses of some superparticles are shown for Mmess = 10^^ and 
10^^ GeV, respectively. One finds that the lightest superparticle among the MSSM particles 

#5ln Fig. |l|, we require that all the soft SUSY breaking mass squared parameters for sleptons be positive 
to make the figure being independent of -Fb. Notice that such a constraint is slightly different from the 
actual experimental constraint based on the mass eigenvalues. (See Figs.|^ and ^.) Constraint on A/,„oss 
is, however, almost unchanged. 
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Figure 2: Mass spectrum of ri, Xi^ xf and ti, with tan/5 = 30, M„ 
M2 = 600 GeV. 



10^2 GeV and 



800 
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Figure 3: Same as Fig. |, except M^ess = 10^^ GeV and M2 = 500 GeV. 
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is almost always a stau, the superpartner of the tau lepton. Note that in the model 1, the 

stau is not stable since the axino is the lightest superparticle and the stau decays to the 

axino. We will discuss phenomenological and cosmological implications to the stau decay. 

In addition, it is also notable that the slepton masses are more enhanced relative to the 

gauginos if we adopt larger value of Mmess- 

In Table ^ we show the superparticle mass spectrum for some typical cases.^ As we 

mentioned, the lightest superparticle in the MSSM sector is the stau. The spectrum is 

rather compact, and squarks and gluinos are much lighter for the same stau mass than in 

the case of gauge mediation [Q. Another interesting point is that the lightest neutralino is 

higgsino-like, which is not achieved in many scenarios of supersymmetry breaking, including 

IT? 



the minimal supergravity model with all the SUSY particles being lighter than ~ 1 TeV 
In fact, one sees that the lightest and the second lightest neutralinos and the lighter chargino 
are quite degenerate in mass, indicating that they are higgsino-like. 

In the minimal anomaly mediation, the sign of the gluino mass is opposite to that of 
the Wino mass. This causes a potential conflict between the muon anomalous magnetic 
moment and Br(6 — * 57) |]TB|]. The result reported by the Brookhaven E821 experiment 
jTSl suggests a deviation from the standard model prediction, indicating a new contribution 
from new physics. If this is real, the new physics should give a positive contribution to a^, 
which constrains the sign of ^hM2 (> 0) in the context of SUSY. This choice also fixes 
the sign of HhM^ (< in the minimal anomaly mediation). This sign plays an important 
role in the SUSY contributions to b ^ 57. With this choice, the charged Higgs loop 
and the chargino/stop loop give the contributions with the same sign, conflicting with the 
experimental constraints. One of the advantages of our scenario is that the gluino mass 
has the same sign as the Wino mass as far as > 3, and thus the two constraints can be 
simultaneously satisfied. 

Figs. ^ and | depict the SUSY contribution to in the model 1 for Mmess = 10^^ and 
10^^ GeV. (For the formula for the SUSY contribution to a^, see |]19|, and for recent works 
see [0.) In general, for a fixed value of of the gaugino mass, the stau mass decreases as 
tan/5 increases. Furthermore, as mentioned before, slepton mass relative to the gaugino 
mass becomes smaller as we adopt smaller value of the messenger scale. As a result, for 
relatively small M^ess; we obtain an upper bound on tan/3 to evade the experimental 
upper bound on the stau mass, unless the overall scale of the SUSY breaking is large 
enough. (End points of the lines in Fig. ^ correspond to such points.) One consequence 
of this fact is that, when the messenger scale is low, the SUSY contribution to the muon 
anomalous magnetic moment is constrained to be relatively small, in particular, compared 
to the deviation observed by the E821 experiment. Once the messenger scale is pushed up, 
however, the stau mass can be large enough and one finds that the requisite contribution 
of 0(10-9) is easily realized. Notice that, for M^^ss = 10^^ GeV, fa = (10^^ - 10^^) GeV 

T^^The Higgs masses are computed by using the effective potential at one loop level. Higher order 
corrections will somewhat reduce the h'^ mass. 

"^^^In the focus point scenario of the supergravity model however, the lightest neutralino becomes 

higgsino-hke pTf. 
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is required since the parameter A has to be smaller than 0(0. 1) in our model. (See Eq. 
( p.lTp .) For such a large value of fa, energy density of the axion field becomes larger than 



the critical density if the standard evolution of the Universe is assumed ||2T[]. For a possible 
cosmological scenario in this case, see the next subsection. 

We also studied the superparticle mass spectrum for the model 2. A typical mass 
spectrum in this case is given in Table 2 as well. The lightest superparticle in the MSSM 
sector is again the stau. Because the axino is heavy in the model 2, the stau becomes the 
LSP in the whole theory. Since such a charged LSP is ruled out cosmologically [^], it must 
decay through i?-parity violation. Furthermore the superparticle masses are somewhat 
spread compared to the previous case, which may cause the naturalness problem. Thus we 
conclude that the model 2 is less attractive than the model 1. 



3.2 Other Issues 

Let us return to the model 1, and consider the effects of the stau decay into the axino. 
The axion-multiplet coupling to the matter is suppressed by the decay constant /„. Since 
we are considering the hadronic axion, the axion coupling to leptons vanishes at tree level. 
In fact, it arises at two- loop order, which originates from the anomalous coupling of the 
axion to two photons. Equivalently one can compute this coupling by using the superfield 
techniques. The axion coupling is then contained in the field dependent wave function 
renormalization of the lepton multiplet. The same technique enables us to compute the 
axino-lepton-slepton coupling in a simple matter as in Ref. fl^. Their result is given for 
the right-handed slepton 

C,a = Q^^aW + h.c., (3.1) 



where 

a 



1 2N5iN5 + 33/5) 



aliMracss) - alimz) . (3.2) 



S7r2 11 

It follows from this that the lifetime of the stau is approximately 

r- = Ni ((X)/10i3GeV)^ (200GeV/m^)^ sec. (3.3) 

Thus in collider experiments, the staus do not decay inside detectors, and they are prac- 
tically regarded as stable particles. Thus highly ionized tracks produced by the heavy 
charged particles will be a signature of this scenario P3 . 



Cosmological implications of our model are also interesting since there exist various 
exotic particles in our model, like axion, axino, and saxino, which may affect evolution of 
the Universe. In particular, in the model 1, saxino may significantly affect the cosmology, 
contrary to the conventional scenario. Therefore, let us comment on this point. 

In the model 1, there are two relatively light scalar fields, which are axion and saxino. 
The initial amplitude of these fields are in general displaced from the minimum of the 
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potential and they start to oscillate as the Universe expands. As a result, we should worry 
about cosmological difficulties possibly caused by these coherent oscillations. The effects 
of the axion field have been extensively studied. 

The effects of the saxino is quite different from those of the axion, since in our model the 
saxino is much heavier than the axion. Consequently, the lifetime of the saxino becomes 
much shorter than that of axion, and hence the saxino may decay before the present epoch. 
The saxino a dominantly decays into the axion pair, and the lifetime is calculated as 



m 



X 



647r (X)2 



- . -1/2 

10 sec X c,-^A-M ^l.r ) ( — ) , (3.4) 
^ ^ 100 TeV ' ^ ° ' ' ^ ^ 



where we approximated that the mass of the saxino is comparable to mx- The lifetime is 
much longer than 1 sec as far as A <^ 1 and F$ ^ 100 TeV. As a result, if the coherent 
oscillation exists in the early Universe, it decays after the big-bang nucleosynthesis (BBN) 
starts. Therefore, if the energy density of a is too large at the time of the BBN, it spoils 
the great success of it. Assuming that the initial amplitude of the saxino field as ~ /a, the 
saxino energy density normalized by the entropy density is estimated as 

P. - -^fl - 1 X 10» GeV X A='= (3.5) 



s mI''^ V 100 TeV 

Notice that this ratio is independent of time. 

Important constraint on the saxino abundance is from the overproduction of ^He. If 
the saxino energy density is too large at the time of the neutron decoupling, it boosts 
up the expansion rate of the Universe and increases the freeze-out temperature of the 
neutron. If this happens, the neutron number density after the freeze out becomes larger, 
resulting in an overproduction of ^He. To avoid this problem, the ratio p„/s should be 
much smaller than the ratio Prad/s at the time of the neutron decoupling, and we obtain 
an upper bound on A of ~ 0(10^"^), corresponding to Mmess ~ 10^ GeV. As seen in Fig. |I|, 
with the messenger scale lower than ~ 10^ GeV, a large number of the messenger multiplet 
should be introduced in order to avoid the tachyonic sfermion mass, which makes the gauge 
coupling constants non-perturbative below the GUT scale. One might think that we may 
have a consistent scenario of cosmology if Mmess ~ 10^ GeV. In this case, however, the 
lifetime of the saxino is as long as ~ 10^ sec and hence the saxino decays after dominating 
the Universe. Since the saxino dominantly decays into the axions, energy density of the 
axion becomes 2 — 3 orders of magnitude larger than that of the radiation and hence there 
exists large extra energy density in the form of the relativistic matter (i.e., the axion). 
First, this fact changes the epoch of the radiation-matter equality to z ~ 10. In addition, 
the axion is a very weakly interacting particle and hence the cosmic density fluctuation for 
the scale which enters the horizon before the equality is washed out by the effect of the free- 
streaming. These effects cause a serious difficulty in the galaxy formation. Furthermore, 
some part of the saxino will decay into the pion pair through the QCD anomaly. Such 
pions dissociate the light elements produced by the BBN and spoil the great success of the 
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standard BBN. Thus, we should conclude that the scenario with Mmess ~ 10^ GeV does 
not work. Therefore, if we assume the conventional scenario of cosmology, our model has 
a serious cosmological disaster. 

This problem is, however, easily solved if a large amount of entropy is produced after 
the onset of the saxino oscillation. With a large entropy production, the energy density of 
the saxino field is diluted and the ratio Pa/ s becomes much smaller than the value given in 
Eq. ( |3.5|) . In the anomaly- mediated models, the moduli fields may be heavy and they can 
become a natural source of the late-time entropy production. It is also important to point 
out that, if a large late-time entropy production occurs, axinos are also diluted. Thus, 
the axino relic density becomes much smaller than the present critical density even if the 
axinos are thermally produced in the early Universe. In addition, if the entropy production 
occurs after the QCD phase transition, coherent oscillation of the axion field is also diluted 



and the cosmological bound on the axion scale /„ can be relaxed In this case, fa much 
larger than 10^^ GeV becomes possible. This fact may be of some help if a relatively large 
value of the SUSY contribution to the muon magnetic moment is required, as mentioned 
in the previous subsection. It should be also noted that the relic staus are also diluted by 
the entropy production. Thus, the relic staus do not significantly affect the BBN even if 
their lifetime is longer than ~ 1 sec. The scenario with the late-time entropy production by 
the modulus decay may be tested by precisely measuring the cosmic microwave anisotropy 
since the decay of the modulus field may generate correlated mixture of the adiabatic and 
iso curvature fluctuations p5|] . 



In the model 2, the axino and saxino are as heavy as ~ 0(-F$). As a result, they are 
short-lived, and irrelevant in cosmology. 



4 Summary 

In this paper, we have considered hadronic axion models in the supersymmetric theory 
where SUSY breaking is mediated by the super- Weyl anomaly. SUSY breaking plays an 
important role in determining the VEV of the axion multiplet X which spontaneously 
breaks f/(l)pQ symmetry. We constructed three models in each of which the role of the 
SUSY breaking is different. The first model has no superpotential for X and thus the scalar 
potential solely comes from SUSY breaking. In fact, the VEV of X is determined by the 
balance between the two terms: one is the negative soft SUSY breaking mass generated 
by the super- Weyl anomaly, and the other is from the higher order term in the Kahler 
potential. In the second model, the PQ symmetry would not be broken in the absence 
of the SUSY breaking. Once the SUSY breaking is switched on, the potential minimum 
shift to the vacuum with the symmetry breaking. In the third model, the PQ symmetry is 
already broken in the exact SUSY limit, but the vacuum is degenerate along a non-compact 
flat direction. It is then the effect of the SUSY breaking that lifts the fiat direction, thus 
fixing the PQ scale. In all three models, the PQ scale can be adjusted to the allowed range 
by appropriate choice of the coupling constants. 
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In the first two models, the PQ multiplet possesses non-trivial SUSY breaking, and 
thus it mediates SUSY breaking to the MSSM sector via a mechanism similar to the 
gauge mediation. Then soft SUSY breaking masses are deflected from those from the pure 
anomaly mediation. Our models provide a realization of the deflected anomaly mediation 



advocated in Ref. On the other hand, in the model 3 the X field has trivial SUSY 
breaking, and thus it does not affect the mass spectrum of the superparticles in the MSSM. 
Therefore the sleptons remain tachyonic unless some other mechanism lifts up the slepton 
masses. 

Phenomenological issues of the models 1 and 2 are briefly discussed. In these models, 
the stau becomes the lightest superparticle in the MSSM sector. In colhder experiments, 
charged tracks due to the stau will be an important signature of the models. We also com- 
puted the SUSY contribution to the muon anomalous magnetic moment. We found that, 
despite relatively heavy Winos, the contribution from the SUSY sector can be comparable 
to or even larger than that from the weak interaction sector. Thus it can easily explain the 
possible discrepancy from the standard model prediction reported by the E821 experiment. 

One of the crucial differences between the two models is the mass of the axino. In the 
model 1, the axino is light and becomes the lightest superparticle of the whole theory. In 
this case, the stau is no longer stable, but decays to the axino, with lifetime typically larger 
than 1 sec. In the model 2, the axino acquires mass at the tree level and thus it is as heavy 
as the gravitino. In this case, the stau will be the LSP of the whole sector and so we need 
i?-parity violation to avoid the charged stable particle. 

We also discussed cosmology of the models. In model 1, the saxino is light and long- 
lived. We argued that its coherent oscillation would spoil the success of the BBN. To 
avoid it, we invoke the late time entropy production, and thus cosmological evolution of 
the model 1 should differ from the standard thermal history of the Universe. On the other 
hand, in model 2, the saxino has a mass comparable to F$, and thus is heavy. Thus it is 
short-lived, and cosmologically harmless. 



Acknowledgment 

The authors thank T. Takahashi for useful discussion. This work was supported in part by 
the Grant-in-aid from the Ministry of Education, Culture, Sports, Science and Technology, 
Japan, priority area (#707) "Supersymmetry and unified theory of elementary particles," 
and in part by the Grants-in-aid No. 11640246 and No. 12047201. 

References 

[1] V. Baluni, Phys. Rev. D19 (1979) 2227. 

[2] R.D. Peccei and H.R. Quinn, Phys. Rev. Lett. 38 (1977) 1440; Phys. Rev. D16 (1977) 
1791. 



18 



[3] S. Weinberg, Phys. Rev. Lett. 40 (1978) 223; F. Wilczek, Phys. Rev. Lett. 40 (1978) 
279. 

[4] J.E. Kim, Phys. Rev. Lett. 43 (1979) 103; M.A. Shifman, A.I. Vainshtein and V.I. 
Zakharov, Nud. Phys. B166 (1980) 493. 

[5] M. Dine, W. Fischler and M. Srednicki, Phys. Lett. B104 (1981) 199; A.P. Zhitnitskii, 
Sov. J. Nud. Phys. 31 (1980) 260. 

[6] K. Rajagopal, M.S. Turner and F. Wilczek, Nud. Phys. B358 (1991) 447; T. Goto and 
M. Yamaguchi, Phys. Lett. B276 (1992) 103; E.J. Chun, J.E. Kim and H.P. Nilles, 
Phys. Lett. B287 (1992) 123; E.J. Chun and A. Lukas, Phys. Lett. B357 (1995) 43; 
T. Asaka and M. Yamaguchi, Phys. Lett. B437 (1998) 51; Phys. Rev. D59 (1999) 
125003. 

[7] L. Randall and R. Sundrum, NucL Phys. B557 (1999) 79; G.F. Giudice, M.A. Luty, 
H. Murayama and R. Rattazzi, JHEP 9812 (1998) 027; J.A. Bagger, T. Moroi and 
E. Poppitz, JHEP 0004 (2000) 009. 

[8] T. Gherghetta, G.F. Giudice and J.D. Wells, NucL Phys. B559 (1999) 27. 

[9] T. Moroi and L. Randall, NucL Phys. B570 (2000) 455. 
[10] A. Pomarol and R. Rattazzi, JHEP 9905 (1999) 013. 
[11] E. Witten, Phys. Lett. B105 (1981) 267. 

[12] L. Randall and R. Sundrum, in Ref. []^; M.A. Luty and R. Sundrum; Phys. Rev. D62 
(2000) 035008. 

[13] See, for example, J. Wess and J. Bagger, Supersymmetry and supergravity (Princeton 
University Press, 1992). 

[14] R. Rattazzi, A. Strumia and J.D. Wells, NucL Phys. B576 (2000) 3. 

[15] J.L. Feng and T. Moroi, Phys. Rev. D61 (2000) 095004. 

[16] J.L. Feng, K.T. Matchev and T. Moroi, Phys. Rev. Lett. 84 (2000) 2322. 

[17] J.L. Feng, K.T. Matchev and F. Wilczek, Phys. Lett. B482 (2000) 388. 

[18] H.N. Brown et al. (Muon g-2 Collaboration), Phys. Rev. Lett. 86 (2001) 2227. 

[19] J.L. Lopez, D.V. Nanopoulos and X. Wang, Phys. Rev. D49 (1994) 366; U. Chat- 
topadhyay and P. Nath, Phys. Rev. D53 (1996) 1648; T. Moroi, Phys. Rev. D53 
(1996) 6565; Erratum-ibid, D56 (1997) 4424. 



19 



[20] J.L. Feng and K.T. Matchev, Phys. Rev. Lett. 86 (2001) 3480; L.L. Everett, G.L. 
Kane, S. Rigolin and L.-T. Wang, Phys. Rev. Lett. 86 (2001) 3484; E.A. Baltz and 
P. Gondolo, Phys. Rev. Lett. 86 (2001) 5004; U. Chattopadhyay and P. Nath, Phys. 
Rev. Lett. 86 (2001) 5854; S. Komine, T. Moroi and M. Yamaguchi, Phys. Lett. B506 
(2001) 93; Phys. Lett. B507 (2001) 224; J. Hisano and K. Tobe, Phys. Lett. B510 
(2001) 197; J.R. Ellis. D.V. Nanopoulos and K.A. Olive, Phys. Lett. B508 (2001) 65; 
R. Arnowitt, B. Dutta, B. Hu and Y. Santoso, Phys. Lett. B505 (2001) 177; K. Choi, 
K. Hwang, S.K. Kang, K.Y. Lee and W.Y. Song, Phys. Rev. D64 (2001) 055001; 
S.P. Martin and J.D. Wells, Phys. Rev. D64 (2001) 035003; H. Baer, C. Balazs, J. 
Ferrandis and X. Tata, Phys. Rev. D64 (2001) 035004; S. Baek, T. Goto, Y. Okada and 
K. Okumura, Phys. Rev. D64 (2001) 095001; C.H. Chen and C.Q. Ceng, Phys. Lett. 
B511 (2001) 77; Z. Chacko and G.D. Kribs, Phys. Rev. D64 (2001) 075015; T. Blazek 
and S.F. King, Phys. Lett. B518 (2001) 109; G.-C. Cho and K. Hagiwara, Phys. Lett. 
B514 (2001) 123; M. Byrne, C. Kolda and J.E. Lennon, |hep-ph/0l08T^ ; S. Komine 



and M. Yamaguchi, |hep-ph/011003^ ; M. Endo and T. Moroi, |hep-ph/0110"383 



[21] J. Preskill, M. Wise and F. Wilczek, Phys. Lett. B180 (1983) 127; L. Abbott and 
P. Sikivie, Phys. Lett. B180 (1983) 133; M. Dine and W. Fischer, Phys. Lett. B180 
(1983) 137. 

[22] A. Kudo and M. Yamaguchi, Phys. Lett. B516 (2001) 151. 

[23] M. Drees and X. Tata, Phys. Lett. B252 (1990) 695; J.L. Feng and T. Moroi, Phys. 
Rev. D58 (1998) 035001. 

[24] M. Kawasaki, T. Moroi and T. Yanagida, Phys. Lett. B383 (1996) 313. 

[25] T. Moroi and T. Takahashi, |hep-ph/0110096| . 



20 



